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A simple, approximate analytical formula is proposed for the reflection function of a semi-infinite, homoge- 
neous particulate layer. It is assumed that the zenith angle of the viewing direction is equal to zero (thus 
corresponding to the case of nadir observations), whereas the light incidence direction is arbitrary. The for- 
mula yields accurate results for incidence—zenith angles less than approximately 85° and can be useful in ana- 
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lyzing satellite nadir observations of optically thick clouds. © 2002 Optical Society of America 


OCIS codes: 010.1290, 280.1310, 290.7050. 


1. INTRODUCTION 


The reflection function R(@), 0, ¢ ) is defined as the ratio 
of the reflected diffused light intensity I( 0), 6, ¢) for the 
case of an arbitrary light-scattering, plane-parallel me- 
dium to that of a Lambertian ideally white reflector. 
That is, we have 
RG 6 I ( bo ? 0, h ) 
( Oa Vy $ ) T® ( bo) , 
where 
I*(09) = F cos 69 
is the intensity of light reflected from the ideally white 
Lambertian reflector, mF is the incident light flux on the 
area perpendicular to the direction of incidence, 0 is the 
light incidence angle, 6 is the observation angle, and ¢ is 
the relative azimuth between incidence and observation 
directions. We also will use 7 = |cos 4, 1 = cos 6. 
Clearly, it follows for the Lambertian ideally white reflec- 
tor thatR = 1. By definition, this result does not depend 
on the viewing geometry. 

Although nonabsorbing optically thick media are white 
when looking from the illumination direction, their reflec- 
tion functions generally differ from 1. The value depends 
on the viewing geometry, the microstructure and thick- 
ness of the medium, and the type of underlying surface. 

We will assume in this study that the optical thickness 
ris such that the reflection function does not change if we 
further increase 7 (i.e., increase the concentration of par- 
ticles or geometrical thickness of a turbid medium). Ef- 
fectively it means that we consider the case of a semi- 
infinite layer. 

The reflection function of a homogeneous, semi-infinite, 
turbid medium plays an important role both in the radia- 
tive transfer theory!” and in planetary optics.® In par- 
ticular, the reflection function of an optically thick scat- 
tering layer can be easily calculated if the reflection 
function of a semi-infinite layer with the same local opti- 
cal characteristics is known.’~° 
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The primary task of this paper is to obtain an approxi- 
mate analytical solution for the reflection function of a 
semi-infinite, plane-parallel layer of cloud water droplets. 
However, the results obtained can also be applied to other 
types of turbid media (especially for media with values of 
the average cosine of the scattering angle that are close to 
those of cloudy media). The reflection function depends 
on the angular scattering coefficient of an elementary vol- 
ume of the turbid medium in question, the viewing and 
observation geometries, and the probability of photon ab- 
sorption in the single-scattering event y. It is defined as 
the ratio C bs/Text, Where oh, and Cext are absorption and 
extinction coefficients of the turbid medium. We will as- 
sume that y = 0, which simplifies the analysis. How- 
ever, the general approach used here can also be applied, 
with a slight modification, to the case of weakly absorbing 
media. Indeed, it is known that)? 


R(n, & $) = Rm, & $) -yK MKA D 


as y— 0. Here R(n, €, ¢) is the reflection function of 
a semi-infinite, weakly absorbing medium, R?(ņ, é, ¢) 
the reflection function of a semi-infinite, nonabsorbing 
medium, y = 4k/3(1 — g), g is the asymmetry param- 
eter, and k is the diffusion exponent of the radiative 
transfer theory (k = [3y(1 — g)]"? as y— 0). It fol- 
lows that Ko(7) = (3/7)(1 + 27) at 7 > 0.2 with the er- 
ror less than 2%.+° Equation (1) allows us to find 
R(n, é @) for a weakly absorbing layer if the function 
R2(n, £, $) is known. 


2. REFLECTION FUNCTION 


The reflection function R°(7, é, &) can be represented as 
the sum of two terms: 


Ram é, $) = Rela £ p) + Radlm& $). (2) 


The function R,.(7, €, @) is the single-scattering contri- 
bution given byf 
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p(y) 
R.A, & 6) = i (3) 


(n+ £? 
where p(®) is the phase function and 
Ə = arccos[ -né + (1 — A1 — &) coso] (4) 


is the scattering angle. 

The function R,,.(7, é& 6) gives us the multiple- 
scattering contribution. Because of the randomizing ef- 
fect of multiple scattering, its dependence on the relative 
azimuth angle and the phase function is weak; so we will 
neglect this dependence. We will use the following ap- 
proximation for this function: 


a + bné+celn + €) 6) 
Rus(y, £) = . 5 
4(7 + &) 
Equation (5) was proposed by Sobolev! for the description 
of the azimuthally averaged reflection function of a semi- 
infinite, nonabsorbing medium. It follows from Eqs. (3)- 
(5) that 


POO) at bné + e(n + é) 
(7 + &) 4(7 + £) 


Rom é ¢) = Z 
(6) 


In order to find the constants a, b, and c, we calculate the 
function R?(7, £, ¢) using the exact radiative transfer 
code.’ The code is based on the numerical solution of 
Ambarzumian’s nonlinear integral equation.’ This al- 
lows us to find the function 


D(m, €) = 4(7 + ORY, £, &) — p8). (7) 
Note that it follows from Eqs. (6) and (7) that 


D(y, €E) =a + byn&+ c(y + €). (8) 


The function D(1, £), calculated with Eq. (7) at 7 = 1 
(nadir observation), is depicted in Fig. 1. Calculations 
were performed for water droplets with an effective 
radius of 6 wm at a wavelength of à = 0.65 um. The 
particle size distribution is given by fr) 
= Nr® exp(—9r/reg), where reg is the effective radius 


@ exact data 
linear fit 
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cosine of the incidence angle 


Fig. 1. Dependence D(1é) as obtained from Eq. (7) (symbols) 
and Eq. (9) (solid curve) for A = 1.468 and B = 7.748. 
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equal to the ratio of the third and second moments of the 
particle size distribution and N is the normalization con- 
stant. The curve in Fig. 1 can be approximated by the 
following linear relationship, 


D(1,) =A + Bë, (9) 


where A = 1.468 and B = 7.748. On the other hand, it 
follows from Eq. (8) for 7 = 1 that 


D1, =a +c + (c +b)é. (10) 


Considering Eqs. (9) and (10), we see that Eq. (10) has 
three free parameters as compared with the two param- 
eters A and B in Eq. (9). We will assume that c = 0. 
Then it follows that a =A and b = B. Equation (6) 
takes the following final form for 7 = 1: 


; a + Bcos bo p(T — 0o) 
Rz( 09) = $ (11) 
1 + cos 6 4(1 + cos 6) 
where 
a = a/4 ~ 0.87, B = 6/4 ~ 1.94. 


Equation (11) reduces the calculation of the reflection 
function of a semi-infinite, nonabsorbing medium at nadir 
observation to the calculation of the phase function p(7 
— o). It should be pointed out that the last term in Eq. 
(11) is small compared with the first term for most inci- 
dence angles. This reduces Eq. (11) to 


0.37 + 1.94 cos A 
R20) = ; (12) 
1 + cos 0o 


which can be used for rapid estimations. 


The comparison of calculations of the reflection func- 
tion for cloudy media with effective droplet radii re = 6 
and 16 um based on Eq. (11) and the exact code described 
in Ref. 7 is presented in Fig. 2(a). One can see that the 
simple approximation given by Eq. (11) can indeed be ap- 
plied to calculations of the reflection function of semi- 
infinite water clouds corresponding to the case of nadir 
observations. The error is less than 2% at 6) < 85° [see 
Fig. 2(b)]. For larger solar angles the accuracy of the ap- 
proximation breaks down as does the assumption of 
plane-parallel geometry. We note also that the reflection 
function of semi-infinite clouds at nadir observation is al- 
most independent of the cloud droplet size, which is not 
surprising given the large size of nonabsorbing cloud 
droplets relative to the wavelength. 

Thus we have managed to parameterize the reflection 
function of a semi-infinite water cloud at nadir observa- 
tions with the simple approximate formula of Eq. (11). 
This allows one to simplify the solution of both the direct 
and the inverse problems of cloud optics, which are based 
largely on the asymptotic theory that is valid for optically 
thick light-scattering layers. For instance, it is known 
that the reflection function of optically thick clouds in the 
visible is given by% 


R(6, 0, $, T) = R2(0, 0, $) — tKo(cos 0)Ko(cos 09), 
(13) 


where 
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—— 16 micrometers, exact data 

F ---O--- 16 micrometers, approximation 
--A--— 6 micrometers, exact data 

: | ---7-- 6 micrometers, approximation 


reflection function 
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Fig. 2. (a) Dependence of the reflection function R(1, £) at nadir 
observation on the incidence angle 4) at \ = 0.65 um and effec- 
tive radii of 6 and 16 wm. (b) Dependence of the error of the ap- 
proximation of Eq. (11) on the incidence angle at à = 0.65 um 
and effective radii of 6 and 16 um. 


1 
t= 3 (14) 
0.75711 -—g) +a 


a~1.07, g is the asymmetry parameter, and 
R°(6, 0), $) is given by Eq. (11). The value of t is the 
global transmittance, which is equal to 1 — r owing to the 
energy conservation law. Here 


2 2T T T 
r= <f asf dé cos of dôg cos OoR(, 0, $, 7) 
T J0 0 0 
(15) 


is the spherical albedo. We see that integral (15), which 
is of paramount importance in climate studies, can be eas- 
ily retrieved from Eq. (13), even if information on the 
cloud optical thickness and microstructure is not avail- 
able. That is, it follows that 
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R?(9, bo, $) = R(9, 0o, $, T) 
r=1 ; (16) 
Ko(cos 0)Ko(cos 6) 


where R(0, 09, $, 7) is the measured function. 
Note that values of 7 and g can be related to the droplet 
size by means of the following approximate equations’: 


3W Iel 1 a7) 
= 1+ i = 0.88 — —;z, 17 
TT apra) R F 2x2 


where X og = 27T, pis the specific density of water, and 
W is the liquid water path. Equations (13), (14), and (17) 
in combination with Eq. (11) can be applied, e.g., to the 
liquid water path determination if the value of the effec- 
tive radius is known (e.g., from measurements of the re- 
flection function in the infrared region of the electromag- 
netic spectrum). 


3. CONCLUSIONS 


The main result of this paper is the simple Eq. (11) for the 
reflection function of a cloudy layer at the nadir observa- 
tion. It can be applied to the remote sensing of cloud mi- 
crostructure from space-based platforms. The weak de- 
pendence of the reflection function described by Eq. (11) 
on the cloud droplet sizes means that one can use this 
function for vertically inhomogeneous clouds without any 
changes. For rapid estimations one can use the simpli- 
fied Eq. (12), which depends only on the geometry of the 
problem. 

I believe also that the results and the general param- 
eterization approach can be applied to turbid media other 
than clouds (especially if they have values of the average 
scattering angles in the single-scattering event close to 
those characteristic for cloudy media). Furthermore, the 
results presented here can be used for the derivation of 
cloud total reflectance from the reflection function mea- 
surements in the visible, which is of importance for cli- 
mate studies. 
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